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Abstract 

We determine a general link between two different solutions of the MaxEnt vari- 
ational problem, namely, the ones that correspond to using either Shannon's or 
Tsallis' entropies in the concomitant variational problem. It is shown that the two 
variations lead to equivalent solutions that take different appearances but contain 
the same information. These solutions are linked by our transformation. 
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1 Introduction 



Nonextensive statistical mechanics (NEXT) [T|2|3] , a generalization of the or- 
thodox Boltzmann-Gibbs (BG) one, is actively investigated and applied in 
many areas of scientific endeavor. NEXT is based on a nonadditive (though 
extensive [3]) entropic information measure, that is characterized by a real 
index q (with q = 1 recovering the standard BG entropy). It has been used 
with regards to variegated systems such as cold atoms in dissipative opti- 
cal lattices [5], dusty plasmas [6], trapped ions [TJ, spinglasses jS], turbu- 
lence in the heliosheath [§], self-organized criticality [TO], high-energy experi- 
ments at LHC/CMS/CERN [IT] and RHIC/PHENIX/Brookhaven [12], low- 
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dimensional dissipative maps [13], finance |14| . galaxies [15], Fokker-Planck 
equation's applications [16j, etc. 



A typical NEXT feature is that it can be expressed by recourse to general- 
izations a la q of standard mathematical concepts [T7]. Included are, for in- 
stance, the logarithm [q-logarithm] and exponential functions (usually denoted 
aS €-q\X ), with e q= i(x) = e x ), addition and multiplication, Fourier transform 
(FT) and the Central Limit Theorem (CLT) [18] . The q- Fourier transform F q 
exhibits the nice property of transforming q-Gaussians into q-Gaussians |18] . 
Recently, plane waves, and the representation of the Dirac delta into plane 
waves have been also generalized [T9f2U|l21|l2"2"] . 

Our central interest here resides in the q-exponential function, regarded as 
the MaxEnt variational solution [3] if the pertinent information measure is 
Tsallis' one. We will show that there is a transform procedure that converts any 
Shannon- MaxEnt solution [23|[2l] into a q-exponential, without modification 
of the associated Lagrange multipliers, that carry with them all the physics of 
the problem at hand. Why? Because of the Legendre transform properties of 



the MaxEnt solutions [see for instance [2^1l2^1l2^1l2Wril2^1l2^^ and references 
therein] . 



Accordingly, we are here proving that the physics of a given problem can be 
discussed in equivalent fashion by recourse to either Shannon's measure or 
Tsallis' one, indistinctly. 



2 The central idea 

We wish to connect orthodox exponentials with q-exponentials. Let us consider 
the Shannon-MaxEnt solution for a constraint given by the average value of 
the variable in question, that we call u: 

p(u)du = exp [—jjL — Xu] du = e _At e~ Xu du, 



where the two Lagrange multipliers /i and A correspond, respectively, to nor- 
malization and conservation of the u— mean value < u >= K. 

Consider now a second variable x such that dx/du = g(x), with g(x) a function 
we will wish to determine below. Assume that in the second variable we express 
the Tsallis-MaxEnt solution, with the same constraints, but employing the 
above mentioned q-exponential functions (e q (x) = [1 + (1 — q)x] 1 ^ 1 ~ q \ The 




(1) 
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support of this function is sometimes finite, depending on the q— value. See 
more details in, for instance, [2]). 

p(x)dx = C e q (—Xx) dx; C = normalization const. (2) 
We want 

p(x)dx = p(u)du. 

This entails: 

Ce q {-Xx)dx = p{x)dx = e-»e- Xu[x) {dx/g{x)), (3) 

that is, given that du = g(x)~ l dx, 

C g{x) e q (-Xx) = exp [~X J dxg{x)~ 1 ]. (4) 

Remember now that 
de q (x) 



dx 



e q (xy, (5) 



so that, taking the logarithm of Eq. (jl]) we find 

In C + In g(x) + In [e q (—X x)] — —fj, — X Jdxg(x)' 1 . (6) 

Now, we derive w.r.t. x Eq. and have 



9{x) 



which leads to a differential equation for our desired g(x): 

g'(x)-Xe q (-Xxy- 1 g(x) + X = 0. (8) 

Solving this equation we establish the link we are looking for. 

3 The differential equation 

For simplicity we now set 

P(x) = -Xe q {~Xx) q ~ l - Q(x) = -A, (9) 
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and recast (jSJ) as 

g'{x) + P{x)g{x) = Q(x). 

Introduce now the integrating factor / 



I = cxp[J dtP(t)}, 



and multiply (|10p by it 

g'(x)I + P(x)Ig(x) = Q{x)I. 



Note that 
d 



[g(x)I]=g'(x)I + P(x)Ig(x), 



dx 

so that (1TTT) becomes 
d 



dx 



[g(x)I] = Q(x)I. 



Integrating this we have now 
g(x)I = Q(x) J Idx + c. 

Finally, we can formally express our "solution" function g(x) as 

/ \ / Q(x)Idx + c 
= ; • 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



Thus, for the differential equation 

g'(x)-\e q (-\x)^g(x) + X = 0, 

we obtain the solution 



(17) 



#(x) = exp (A J egi-Xx'Y^dx' 

xlx 1 

-X I exp -A J e^-Xx'^-^dx" \ dx' + c 



where c is an integration constant. Now, using 

X 

X J egi-Xx'^dx' = In (e q (-Xx) { ~^ 



and 



/ x' 



— X exp 



A / e q (-Xx")^dx" dx 



e q (-Xx)^ 
2-q 



we finally arrive at 



g{x) = e q {-Xx)^ 



e q (-Xx)^ 



2-q 



+ c 



which provides us with the Tsallis-Shannon Jacobian 
1 



J(x) 



9{x) 



Consider now the instance q — > 1. One obviously ought to have g(x 
face 



g{x)^e 



Xx 



e~ Xx + c 



^Xx 



= l + ce" 
which entails c = and 



9K X ) 



e q (-XxY 1 ^ l-(l-q)X 



x 



2-q 



2-q 



3.1 Expansion near q=l 



Let us now take q = 1 — e in the solution g(x) = g(x;q) (Eq.(l25 

1-eXx 

g{x;q = 1 - e) 



1 + e 



Thus, a first-order expansion in e gives: 



g(x; q = 1 - e) = 1 - (1 + Ax) e + 0(e 2 ). 



(27) 



3.2 TTie case g = 2 



Now, let us look at the q = 2 — e scenario for our solution g(x) = g(x; g); 

l + (l-e)Ax 

2(x;g = 2-e) = ^ '—. {21 



That, in the limit of e — > we have 



(1 + Ax) , , 

s(x;g = 2-e)x^ 29 

e 

->+oo as e -> + (i.e.,g->2~) and (30) 

->-oo as e -> 0~ (i.e., g -> 2 + ). (31) 

There is a divergence in the first term, in the form of 1/e. 

We are then in a position to state that, for Ax > —1, and given the definition 
of the q-exponential (see [2]), our "inverse- Jacobian" function g(x) is positive 
for q < 2 and negative for q > 2, while diverging at q — 2. We conclude that 
the transform we are studying is not valid only in the isolated case q = 2 and 
changes sign there. 



4 Arbitrary constraint 



We generalize now the preceding considerations to the case of a generalized 
constraint < h(x) >, with h 6 £ 2 - The concomitant Shannon MaxEnt solution 
is 



p{u)du = e-^e- Xh(u) du, (32) 

(h(u)) = K, (33) 

p{u)du = l. (34) 
while the Tsallis-MaxEnt solution with the same constraint becomes 

p(x)dx = C e g (— A h(x)) dx; C = normalization const. (35) 
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Assume that u{x) exists and call, as before, dx/du = g(x). We have, as our 
cornerstone the relation 



p(x)dx = p(u)du, 



(36) 



entailing 

C e g (-X h(x)) dx = e - " e~ Xh{u{x) \dx/ g{x)), 



(37) 



so that, taking the logarithm to this equation we find 
In C + In [g(x)} + In [e q (-X h(x))} = -fj, - X h(u(x)). 



(38) 



Taking derivatives w.r.t. x yields 

44 - A eJ-X hix^h'ix) = -X -^h'(x), 
g(x) g{x) 



(39) 



which leads to a differential equation for our desired transformation function 



g\x) - X e q (-X h{x)) q ^ h\x)g{x) + Xh'(x) = 0, 



(40) 



quite similar in shape to Eq. ffTTl) . being thus solved in similar fashion. Using 
Eq. ( [161) we encounter 



P = j e q (-Xh(x)) q ^(-Xh\x))dx, 

= ln(e 9 (-A/i(x))), 

e-f p = e g (-Xh(x))-\ 

Qel p = J e q (-Xh(x))(-Xh'(x))dx, 

_ e q (-Xh(x)) 2 - q 
2-q ' 



that leads to 

g(x) = e g (-A/i(x)) _1 



e q {-Xh{x)) 2 ~i 
2-q 



+ c 



(41) 
(42) 
(43) 
(44) 

(45) 



(46) 



It is clear that for q = 1 we have g(x) = 1 if c=0. Thus, 
1 - (l-q)Xh(x) 



2-q 



(47) 
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Also, the regime-change at q = 2 discussed above does not change. 
4-1 Special case h(x) = x 2 

For the special case of a variance constraint we have 

g'(x) - X e q (-X x 2 ) q ~ l 2xg(x) = -X2x, (48) 

whose solution is 

I — (I — q)Xx 2 , t , 

z q 

Let us take now q = 1 — e in the solution g(x) = g(x; q). Then, an expansion 
near e = yields 

g{x; q = 1 - e) = 1 - (l + Ax 2 ) e + 0(e 2 ). (50) 

5 Generalization to M constraints 

Now, we generalize to the case of M constraints of the form: 

(h l (u)) = K i ; z = l,...,M. (51) 

Lets use vector notation and call h = (hi, h 2 , . . . , Hm) and A = (Ai, A2, • • • , Xm) 
its associated Lagrange multipliers. Then we have for the Shannon MaxEnt 
solution 

p(u)du = e-> 1 e- Xh{u) du, (52) 
while the Tsallis-MaxEnt solution with the same constraints becomes 

p(x)dx = C e q (—X ■ h(x)) dx; C = normalization const. (53) 

Assume that u(x) exists and call, as before, dx/du = g(x). Then, following 
the steps of the previous section and solving the corresponding differential 
equation, we obtain 



g(x) = e q (-X-h(x)r 1 



e q (-X-h(x)) 2 -* 
2-q 



(54) 
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where c = in order to get g(x) — 1 in the limit q — > 1. Thus, 



1 - (1 - q)X ■ h(x) 
2^ 



(55) 



6 Conclusions 

We have shown here that, from a MaxEnt practitioner view-point, one can 
indistinctly employ Shannon's logarithmic entropy S or Tsallis' power-law 
one S q (for any q except q = 2). The physics described is the same. To choose 
between S and S q is just a matter of convenience in the sense of getting simpler 
expressions in one case than in the other. 

The link between the two concomitant probability distributions Pshannon{%) 
and PTsaiHs( x ) is given by the Jacobian J = 1/g, where g is the simple function 



with A the pertinent Lagrange multiplier, and h(x) G £2 an arbitrary function 
whose mean value < h > constitutes MaxEnt 's informational input. 



7 Appendix: The four different Tsallis' treatments 

A savvy Tsallis practitioner may wonder what happens with the four different 
ways of computing q— mean values that one finds in Tsallis' literature (see 
|31] and references therein). In addition to the normal expectation values we 
have employed above, one also encounters, for a quantity A(x) averaging ways 
that, themselves, depend upon q (see below). This transforms our cornerstone- 
equality pshannon{u)du = pTsaiUs{x)dx into something much more complicated. 
However, there is a way out, following the discoveries reported in |31) . 

Bernoulli published in the Ars Conjectandi the first formal attempt to deal 
with probabilities already in 1713 and Laplace further formalized the subject 
in his Theorie analytique des Probabilites of 1820. In the intervening centuries 
Probability Theory (PT) has grown into a rich, powerful, and extremely useful 
branch of Mathematics. Contemporary Physics heavily relies on PT for a large 
part of its basic structure, Statistical Mechanics [32 33 34J, of course, being a 
most conspicuous example. One of PT basic definitions is that of the mean 



1 - (1 -q)Xh(x) 
2^ 



(56) 
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value of an observable A (a measurable quantity). Let A stand for the linear 
operator or dynamical variable associated with A. Then, 



(A) = / dxp(x) A(x). (57) 



This was the averging procedure that Tsallis used in his first, pioneering 1988 
paper [1], and the one discussed in the preceding Sections. It is well known 
that, in some specific cases, it becomes necessary to use "weighted" mean 
values, of the form 

(A)= [ dxf\p(x)]A(x), (58) 



with / an analytical function of p. This happens, for instance, when there is a 
set of states characterized by a distribution with a recognizable maximum and 
a large tail that contains low but finite probabilities. One faces then the need 
of making a pragmatical (usually of experimental origin) decision regarding / 
|31j . In the first stage of NEXT-development, its pioneer practitioners made 
the pragmatic choice of using "weighted" mean values, of rather unfamiliar 
appearance for many physicists. Why? The reasons were of theoretical origin. 
It was at the time believed that, using the familiar linear, unbiased mean 
values, one was unable to get rid of the Lagrange multiplier associated to 
probability-normalization. Since the Tsallis' formalism yields, in the limit q — > 
1, the orthodox Jaynes-Shannon treatment, the natural choice was to construct 
weighted expectation values (EVs) using the index q, 

(A) q = J dxp(x) q A(x), (59) 

the so-called Curado- Tsallis unbiased mean values (MV) j35]. As shown in [3], 
employing the Curado- Tsallis (CT) mean values allowed one to obtain an an- 
alytical expression for the partition function out of the concomitant MaxEnt 
process [3J. This EV choice leads to a non extensive formalism endowed with 
interesting features: i) the above mentioned property of its partition func- 
tion Z, ii) a numerical treatment that is relatively simple, and iii) proper 
results in the limit q — > 1. It has, unfortunately, the drawback of exhibiting 
un-normalized mean values, i.e., ((1)) ^ 1. The latter problem was circum- 
vented in the subsequent work of Tsallis- Mendes-Plastino (TMP) [3H], that 
"normalized" the CT treatment by employing mean values of the form 

(A) q = J dx^^A(x); X q = J dxp(x) q . (60) 

Most NEXT works employ the TMP procedure. However, the concomitant 
treatment is not at all simple. Numerical complications often ensue, which has 
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encouraged the development of different, alternative approach called the OLM 
one |37| . that preserves the main TMP-idea (the X q normalization sum) but 
is numerically simpler. Now, despite appearances, the four Tsallis' treatments 
are equivalent, as shown in [3T]. By equivalence we mean that if one knows 
the probability treatment Pf, i = 1,2,3,4 obtained by anyone of the four 
treatments, there is a unique, automatic way to write down Pj\ j ^ i. More 
precisely, any Pj is a q-exponential, and they all possess the same information- 
amount if the pair q, /3 is appropriately "translated" from a version to the other 
HSU- Indeed, 

Z = 5>xp 9 ,(-/3*x). (61) 

i 

Then, as made explicit in [31| . given any of the four possible Pj's q*, /3*, 
one can get the q, (3 values appropriate for the q-exponential of any j ^ i. 
Such "dictionary" allows one to translate the results obtained in the preceding 
Sections to any other of the three remaining averaging procedures. 
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